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Abstract: The rate of events measured with the surface detector of the Pierre Auger Observatory is found 
to be modulated by the weather conditions. This effect is due to the increasing amount of matter traversed 
by a shower as the ground pressure increases and to the inverse proportionality of the Moliere radius to 
the air density near ground. Air-shower simulations with different realistic profiles of the atmosphere 
support this interpretation of the observed effects. 



Introduction 

The surface detector (SD) of the Auger Southern 
Observatory, located in Malargiie, Argentina, is de- 
signed for the detection of ultra high energy cos- 
mic rays through the measurement of the signal in- 
duced by the shower particles reaching the obser- 
vation level (~ 880 g cm~ 2 ) in an array of water- 
Cherenkov tanks arranged in a triangular grid with 
1500 m spacing. 

The regular data taking of the SD started in January 
2004, with the array continuosly growing from 100 
stations up to the current 1200. To check the de- 
tector stability we monitored the shower counting 
rate, finding that it is modulated by weather effects. 
This dependence is expected since changes in the 
atmospheric density profile due to weather varia- 
tions influence the development of the air shower 
and in turn the amplitude of the signal measured at 
ground. As a consequence, a study of the detector 
stability has to account for the rate dependence on 
the atmospheric conditions. Moreover, since the 
SD estimate of the energy of the primary particle 
is based on 5(1000), the signal measured at 1000 
m from the shower axis, we are interested in the 
dependence of 5(1000) on the atmospheric con- 
ditions. This requires a continuous monitoring of 
the weather and a good knowledge of the relation- 
ship between 5(1000) and the measured weather 
parameters. The former is provided by a meteo- 



rological station, at the centre of the SD array, that 
records the weather parameters every 5 min, allow- 
ing the correlation of the modulation of observed 
quantities, such as the rate of events, with the mea- 
sured ground temperature T and pressure P. 

Weather effects on EAS 

The expected effects related to the change of 
weather conditions are essentially two: 

(i) an increase in the ground pressure P corre- 
sponds to an increased slant depth X and implies 
that the shower is older when it reaches the ground 
level. The longitudinal development of the elec- 
tromagnetic component of the shower at 1 km 
from the core can be parameterised as a Gaisser- 
Hillas profile, N em (E, X) cx X k ™' K exp[(X TO - 
X)/A], where X m is the average maximum of the 
shower at 1 km from the core (~ 200 g cm -2 
deeper than at the core) and A ~ 70 g cm -2 is an 
effective hadronic attenuation length. Then, under 
a pressure change, the electromagnetic component 
S em of 5(1000) changes by 



d ln5 e 
dP 



~X~ 



sec 9 



(i) 



where dX = dP sec 6 was used. Since for the en- 
ergies of interest, E > 10 18 eV, the maximum of 
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Figure 1 : Daily averages of ground pressure (top), 
density (middle) and event rate (bottom, black). 
Since the pressure is stable, the prominent effect 
on the rate modulation is due to the density (tem- 
perature) variation. The red points in the bottom 
plot show the results of the fit. 

vertical showers is close to ground, this effect is 
expected to be more pronounced for inclined show- 
ers. 

(ii) an increase in the air density reduces the 
Moliere radius tm (proportional to 1/p) and hence 
the lateral extent of the electromagnetic compo- 
nent of the shower. The lateral distribution of 
the electromagnetic component can be approx- 
imately described with an NKG profile, which 
for large radius r from the core behaves as 
N em (r) oc rVf (r/rjvf)~ a , where a ~ 4 and 
tm ~ 83 m/(yo/ kg m~ 3 ). Hence, under a den- 
sity change 

d lng em ^ (2 - a) 
dp p 

The effective value of tm is that corresponding to 
the air density p* two cascade units above ground 
[6] (~ 700 m cos 8 at the Auger site, with 8 being 
the zenith angle). Since the ground T, P are the 
only available observables, we have to express p* 
in terms of the density p measured at ground. 



Figure 2: Variation of P (top) and p (middle) dur- 
ing the day: the values averaged over 2005 and 
2006 (lines) are shown together with the maxi- 
mum variation values during the 2 years consid- 
ered (crosses). Bottom: the result of the fit (red) 
reproduces very well the average diurnal modula- 
tion of the measured rate (black). The local time is 
UTC - 3 h 

On time scales of one day or more, the temperature 
gradient in the lowest layers of the atmosphere (the 
planetary boundary layer) can be described by an 
average value of 6.5°C km -1 ; therefore the varia- 
tion of p* is the same as that of p. An additional 
effect is related to the diurnal variations of the gra- 
dient that is smaller before sunrise, at which time 
even T inversions are common, and larger in the 
early afternoon hours. As a result, the amplitude 
of diurnal variations in T (and p) is smaller at 2 cu 
than at ground level by a factor ~ 0.5. We define 
the average daily densities pd and p* d and the ref- 
erence values (averaged over 2 years of measure- 
ments) p = 1.055 kg m~ 3 and P = 861.9 hPa, 
Pq denotes the reference density at 2 cu above 
ground. The energy reconstructed with no cor- 
rection for weather effects is E r oc [S^IOOO)] 5 , 
where B = 1.13 ± 0.02 [1]. Hence we can pa- 
rameterise the relation between the shower energy 
Eq (8, P, p) at the reference weather conditions and 
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the reconstructed one E r a.si 

E = E r {l - a P (P - P ) - a p (p* - p* )} B 
= E r {l - a P {P - P ) - 

a P {pd - Pa) - P P (p - Pd)} B (3) 

where the coefficients a p ,p and f3 p depend on the 
zenith angle 9. 

Assuming that the cosmic ray spectrum is a pure 
power law d J/dE oc E^ 1 ', it is easy to show that 
the rate R(0, P, p) of events at a given zenith angle 
9 can be expressed as: 

R = R {1 + a P {P - P ) + 

a p (pd - Po) + b p (p - p d )} (4) 

with Po = R(9,Po,Po) and coefficients a p .p = 
(Bj - l)a p ,p and b p = (Bj - l)/3 p , the latter 
describing the diurnal modulation of the rate with 
the density. 

Modulation of the measured rates of 
events 

To study the modulation of the event rate with the 
ground weather parameters, we use the data taken 
from 1 January 2005 to 31 December 2006 that 
have a zenith angle 9 < 60°. The data selection 
criterion is the same as applied for the SD spec- 
trum [ 1 ] . The value of the air density p at ground 
is deduced from P and T measured at the central 
meteorological station. Rather than using the raw 
number of triggering events, we compute the rates, 
as a function of time, to account for the temporal 
variation of the active detection area due mainly to 
the deployment of new stations and occasionally to 
stations experiencing a temporary failure [3]. The 
modulation of the rate during the year, and as a 
function of the hour of the day, follows the changes 
in density and pressure (Figs. 1 and 2). A charac- 
teristic of the Malargtie site is the stability of pres- 
sure (less than ±2% variation), while pd changes 
up to a maximum of ±6% during the year with 
an additional diurnal variation of density of ±2% 
on average, with maximum values of +g% during 
the two years considered. Assuming that the rates 
computed each hour follow a Poisson distribution, 
a maximum likelihood fit gives the estimated val- 
ues of the coefficients in eq. (4) averaged over the 



event distribution in the zenith range 0° — 60° : 

a P = (-0.0009 ± 0.0005) hPa" 1 

a p = (-2.68 ± 0.07) kg^ 1 m 3 (5) 

b p = (-0.85 ±0.07) kg" 1 m 3 . 
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Figure 3: Seasonal atmospheric profiles derived 
from the parameterisation of radio soundings per- 
formed in Malargue and used in simulations. The 
density profiles are divided by the profile of an 
isothermal atmosphere (with Xq = 900 g cm -2 ) to 
enhance the differences. The corresponding values 
of P and T are given in the box. 



Comparison of the experimental results 
with model and simulations 

To test the validity of our interpretation, we com- 
pare the coefficients obtained from the fit of data 
with results from full shower simulations and the 
predictions on theoretical grounds. 
The Corsika code [5] with the QGSjetll model [7] 
for high energy hadronic interactions, was used to 
simulate a set of proton showers at 10 19 eV in 5 
different atmospheres and at various zenith angles. 
The atmospheric profiles used (Fig. 3) are a param- 
eterisation of the seasonal averages of several ra- 
dio soundings carried out at the detector site [4] 
and provide a sample of realistic conditions above 
the Auger SD array, but, being averages on large 
time scales, do not account for the diurnal varia- 
tion of the temperature in the lower atmosphere. 
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The expected signal 5(1000) is estimated through 
the simplified assumptions that e + , and pho- 
tons deposit all their energy in the surface detector, 
while for muons we take the minimum between the 
kinetic energy and 240 MeV (the energy deposited 
by a vertical muon crossing a SD tank). As ex- 
pected, the simulated signal depends on the ground 
density and pressure according to the expression in 
eq. (3) (with p — pd = 0) with coefficients a p and 
up shown in Fig. 4 for all zenith angles between 
0° and 60°. The large uncertainties are due to the 
limited number of atmospheric profiles used. 
For the theoretical expectations, we consider the 
variation of the total signal, given by the sum of the 
electromagnetic and muonic component. The coef- 
ficients a p .p in (3) result from the variation of both 
components: a PtP = F em a e p m P + (1 - F em )a^ P , 
where F em is the electromagnetic fraction at 1 km. 
The dependence of S em on p and P is discussed 
in section 2. For a quantitative prediction we adopt 
in eq. (1) X m = 950 g cm" 2 , typical of 10 EeV 
proton showers, and X = 880 sec 9 g cm -2 . For 
the electromagnetic fraction F em we use a fit to the 
results of shower simulations with 10 EeV protons 
(F em ~ 0.7 near the vertical and decreasing with 
9 to reach ~ 0.2 at 60°). We assume a negligi- 
ble correlation of S p with pressure and a constant 
value at = —0.26 kg -1 m 3 for the dependence on 
density (suggested by the results of simulations). 
In Fig. 4 we compare the coefficients obtained by 
fitting the data in five zenith ranges. The procedure 
to obtain a p p is the same described in section 3, 
then we derive the signal coefficients a p ,p divid- 
ing by (B7 — 1) = 2. Their values are in good 
agreement with both the model predictions and the 
results from simulations. 

Conclusions 

The modulation of the event rates measured by the 
Auger SD can be explained by known effects on 
the shower development, both on seasonal and di- 
urnal scales. At the Auger site the dominant ef- 
fect is related to the density (temperature) vari- 
ation. The systematic error, when determining 
the energy of a single shower in the zenith range 
0° — 60°, amounts to a maximum of ~ 10% (for 
extreme values of ground pressure and tempera- 
ture). The quantitative agreement of the theoreti- 



cal model with simulations and data, suggests that 
it can be used to correct the SD energy reconstruc- 
tion for weather induced effects. 
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Figure 4: Comparison of the signal coefficients a p 
(squares) and a p (circles) obtained from the fit of 
simulated signal for 10 19 eV proton showers, fit 
of the measured rates (shaded rectangles), and val- 
ues obtained with the theoretical model described 
in the text (lines). The value of X m = 750 g cm -2 , 
used in the model, corresponds to 10 EeV showers 
according to the measured elongation rate [2]. 
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